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Abstract: In the inverse seesaw mechanism where the spontaneously broken B-L symme- 
try induces tiny B-L violating Majorana masses for right-handed neutrinos, non-standard 
Higgs signatures can arise due to a possible Higgs doublet and singlet mixing and/or Higgs 
boson decays to a left- and right-handed neutrino. This leads to a remarkable feature of 
hadronically quiet di-lepton final states which can exhibit, in particular, lepton flavour 
violating signatures coming from flavour-dependent neutrino Yukawa couplings. In this 
process, one lepton coming from the right-handed decay could be soft enough can be 
missed by the trigger level cuts of CMS and ATLAS for the di-lepton plus missing energy 
signature. The prospects of such a signature are investigated for 8 TeV and 14 TeV center 
of mass energy of the LHC, taking the maximum value of the allowed neutrino Yukawa cou- 
pling and the right-handed neutrino mass of 100 GeV. A PYTHIA level simulation shows 
that the integrated luminosity of 10-20/fb and 1.6/fb for 8 TeV is required to observe the 
inclusive leptonic and lepton flavour violating signatures, respectively. For 14 TeV, the 
reach is more and a larger parameter space of the inverse seesaw model can be probed. 
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1. Introduction 

A Higgs-like boson around 125 GeV has been observed by both CMS f| and ATLAS @ 
collaborations at the LHC. The 5a discovery reach has been obtained for the Higgs boson 
to 77 mode for both CMS and ATLAS. The Higgs boson, if it is Standard Model (SM) like, 
can decay to di-leptonic final states through WW* decay mode, which is of our interest in 
this paper. However, the significance in such a channel is still low: around 1.6a for CMS 
Q and 2.8a for ATLAS ||. The next question is, if such data will be able to provide 
information of the Higgs boson property, whether it is exactly of the SM or deviates from 
the SM. Non-standard Higgs structure arises in many extensions of the SM motivated by 
various theoretical reasoning and/or experimental requirements. 

In the inverse seesaw mechanism Q explaining tiny neutrino mass, the seesaw sector 
couples to the SM sector with a Yukawa coupling of order-one and thus can lead to observ- 
able signatures at the LHC. A promising way to realize the inverse seesaw mechanism is to 
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consider a gauged B — L symmetry which is broken spontaneously at TeV scale J|, |(| . In 
this case, the B — L Higgs boson can mix with the SM Higgs boson which can change the 
SM Higgs boson signature at the LHC. Although the current LHC data strongly restricts 
such a mixing, we will try to take a conservative limit to maximize the mixing effect. Then 
the Higgs sector of the inverse seesaw model contains two Higgs bosons, a light and heavy 
one denoted by h and H, respectively. As these states are mixtures of the doublet and 
singlet Higgs states, their production is governed by the usual gluon fusion process of the 
doublet state. 

On the other hand, the Higgs bosons can decay to the left- and right- handed neutrinos 
(denoted by v and \E r , respectively): 



or H — > hh depending on their masses. Then, the right-handed neutrino allows the decay 
modes: 



which can lead to hadronically quiet final states from the Higgs boson decay. A careful 
observation shows that the lepton coming from the decay of the right-handed neutrino, 
could really be soft, depending on the mass difference between the right-handed neutrino 
and the W boson. We stress that, for some parameter points, this soft lepton could be 
already missed by ATLAS and CMS due to hard trigger level cuts for the leptons in the 
study of the Higgs boson to WW* Q . Another remarkable feature of the inverse seesaw 
is the possibility of lepton flavour violation as ^ can dominantly decay to a specific lepton 
flavour, which occurs generically with hierarchical neutrino Yukawa couplings. In this 
paper, we analyze such non-standard Higgs signatures to discuss its discovery potential at 
the 8 and 14 TeV LHC. 

The paper is organized as follows. In Section 2, we introduce the inverse seesaw model 
presenting the mass spectrum and the interaction vertices. The explicit form of all the 
interaction vertices of the model is shown in Appendix A. In Section 3, we discuss exper- 
imental constraints on the neutrino Yukawa coupling coming from the lepton universality 
as well as on the Higgs boson masses coming from the LEP and current LHC data. In 
Section 4, the benchmark points for our study are set up to calculate the branching ratios 
for the allowed decay model, and the final state phenomenology of each benchmark point 
is discussed. In Section 5, we analyze the missing and lepton pr distributions and lepton 
multiplicity to suppress the standard di-boson background. Then, we discuss the discovery 

1 For related studies, see the recent works IT], H]. 
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potential of the signal events of two lepton (plus missing px) final states at the 8 and 14 
TeV LHC. We conclude in Section 6. 



2. Inverse seesaw model with XJ(1) b-l 



A natural way to realize the inverse seesaw mechanism is to introduce a B — L gauge 
symmetry whose spontaneous breaking at the TeV scale generates a tiny B — L breaking 
Majorana mass for right-handed neutrinos. The "minimal" field content implementing this 
idea would be as follows: 



Particle 
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u c ,d c 


L 


e c ,N 


Si 


s 2 
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Y B -L 


1/3 


-1/3 
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-1/2 
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(2.1) 



where a pair of fermionic S\ and S2 is required to cancel the anomaly. The SM and B — L 
Higgs bosons are denoted by <& and x, respectively. The gauge invariance of the model 
allows the leptonic sector Lagrangian in terms of the Weyl fermions as follows: 

-C = y e L1>e c + y v L&N + ysNxS, + ^ 2 Sf + ^% 2 S 2 2 + h.c. , (2.2) 

where $ c = e<I>* and A is a cut-off scale. Note that the mass term S1S2 can be suppressed 
by introducing, e.g., a discrete symmetry Z2 under which S2 is odd and the others are even. 
After the symmetry breaking, x = (x° + V ')/V^ and $ = (4> + , 4>) T with <p = (0° + v)/y/2, 
the neutrino sector can be written as 

C = m D u'N + M N NS 1 + fi S Sf + h.c. , (2.3) 

where mo = Mjy = ysv'/y/2 and fis = \s x v' 2 /2k. In the flavour basis {v\ N, S\}, 

the 3x3 neutrino mass matrix of one generation takes the form: 

/ m D \ 

m D M N . (2.4) 
\ M N us J 

From the diagonalization of this mass matrix, one can find the light neutrino mass 



(2.5) 



Considering the constraint on tti^/Mtv, which will be discussed in the following section, 
one can obtain m u « 0.1 eV with the cut-off scale A « 10 15 GeV. 

For the consideration of collider phenomenology, one can simply take the limit of 
Us = 0. In this case, the active neutrino becomes massless and the other two heavy (Weyl) 



-3- 



neutrinos are degenerate to form a Dirac fermion. The mass basis is defined by the following 
rotation 



Si 



cos a sin t 
— sin 9 cos I 



v 

N c 



where sin 6 



m D 



m 2 D + M 2 N 



(2.6) 



where (N, N c ) form a Dirac fermion denoted by \& with the mass = \lm? D + M^. 
Now let us analyze the Higgs potential given by 



V($,x) = m\&$> + 



+ Ai($t$)2 + A 2 ( X f x) 2 + \^^){X ] X)- 



(2.7) 



After spontaneous breaking, the mass matrix is given by 

/ 2X lV 2 X 3 vv' ' 



M(/,^ 



where we use the stationary conditions: 

dV 

dV_ 

dx 



X^vv 1 2A 2 v 



ml + Xiv 2 + ^-X 3 v 



1 



-> m\ + A 2 ?/ 2 + -A3?; 2 



0. 
0. 



(2.9) 
(2.10) 



The neutral component of x an d ^ mixes as follows: 



X 



cos a sin a 
— sin a cos a 



ft 



(2.11) 



where h (H) is the light (heavy) Higgs boson. The mass eigenvalues and the mixing angle 
are given by 



ml jH = Aiw 2 + X 2 v' 2 t V (Aiv 2 - A 2 r/ 2 ) 2 + AlW 2 , 

A3W' 



tan 2q 



A 2 t> /2 — Aif 2 



(2.12) 
(2.13) 



where v = 246 GeV. The most stringent bound on the value of the B—L symmetry breaking 
scale v 1 comes from the LEP II data requiring mz' / 9B-L = \^b—l\ v> ^ ^ TeV [10[. This 
tells us that v' > 12 TeV in our case. We set v' = 12 TeV for our analysis. 

The interacting terms from the Higgs kinetic Lagrangian are given by 



C K = \D^\ 2 + \D lxX ? 



m 2 w W+W^ + -m 2 z Z^ 



1 + 



c a h + s a H 



(2.14) 



1 + 



-s n h + c n H 
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where mz = gv/(2cg w ) and mz 1 = \^b-l\9b-l v ' ■ The leptonic interaction Lagrangian is 
given by 



C 



int. 



= (h cos a + H sin a) q Ui PRq Ui + -j= (h cos a + H sin o^Pr^ 

+ (h cos a + sin a) [cos Q^>iP L Vi + sin 0#jP L *j] 
V2 

+ ^ (/j cos a + H sin a)liP R li + ^j=(-h sin a + # cos a) [- sin 0^>iP L i>i + cos e^iPi*,] 
\/2 V 2 

+ -^W+ [cos0^Pi(^ MJ v5)^ + sin0*n M ^(^MiV5)ii^] 
.9 



- ( _ " ^ [cos 2 OvrfPLVi + sin 2 e^^^Pi^i + cos sin 0(^7^^ + ^ 7 M ^)] 
zcosty^/ y 

+ /i.e. , (2.15) 



where i = (e, fj,, r). and (UMNs)ij is Maki-Nakagawa-Sakata matrix [11|. All the interaction 
terms relevant for our analysis are summarized in Appendix A. 



3. Experimental constraints on y v and the Higgs masses 

The inverse seesaw model allows large neutrino Yukawa couplings y v . The most stringent 



constraint on y u comes from the electroweak precision data. From Table 8 of Ref. [12], one 
gets 



(m D )e iM)T 



M N 



< (0.055, 0.057, 0.079) . 



(3.1) 



This tells us that (y u )i < QMjv/174GeV where q is the value for each flavour given in 
Eq. ( |3.1| ). To be consistent with this bound, we will take 

(m D )i 



smi 



M N 



0.05 



(3.2) 



for i = e or /i for the collider study in the following sections. 



The extra Higgs boson mass is constrained by the LEP and LHC data [13] as well as the 
current LHC data (2|, [l4|]. In the inverse seesaw model, the Higgs production cross sections 
are suppressed by a factor of c 2 and s 2 for h and H, respectively. In addition, the decay 
branching fraction into SM particles is suppressed by a factor of 1 — T$r(h/H — > non-SM). 
If there exists a Higgs boson lighter than the 125 GeV Higgs, it is more singlet-like and 



an upper bound on c 2 is put by the LEP data [13]. For the Higgs boson mass of 50 
GeV, we find c 2 < 0.05. For a Higgs boson heavier than 125 GeV, the current LHC data 
|2|, [14|] puts an upper bound on s 2 as shown in Fig. |l[ For the 125 GeV Higgs boson, 
the overall signal strength measuring the deviation from the SM prediction is found to be 
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150 200 250 300 350 400 450 

rn H (GeV) 

Figure 1: The heavy Higgs boson mass versus s„. The red region is excluded by ATLAS, and the 
blue one is excluded by CMS. The gray region is excluded by the observation of the 125 GeV Higgs 
at CMS. 




Figure 2: Feynman diagrams of the Higgs production via gluon fusion and the decay of Higgses. 

H = 1.4 ± 0.3 (ATLAS) § and \i = 0.87 ± 0.23 (CMS) Q. Thus a non-standard Higgs 
boson contribution to the 125 GeV Higgs is strongly disfavoured. For our analysis, we will 
put a very conservative bound: [i > 0.64, or s\ < 0.36 to see the maximized mixing effect. 
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Figure 3: Standard Model Higgs production via gluon fusion at large Hadron collider for ECM=7, 
8, 14 TeV. CTEQ5L PDF and Q = have been used. For our Higgs boson, we have to multiply 
by the corresponding rescaling factor sin 2 a or cos 2 a. 
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Figure 4: Variation of 'J decay branching fraction with m^. 



4. Benchmark Points and final state phenomenology 



In this article, we are focusing on the gluon fusion process leading to a light/heavy Higgs bo- 
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son and its further decay to the right-handed neutrino and the corresponding phenomenol- 
ogy as stated in the introduction. Fig. || shows the Feynman diagram responsible for this 
process. Fig. ^ describes the variation of the Higgs production cross-section depending on 
the Higgs boson mass for 14 TeV and 8 TeV center of mass energy (ECM) at the LHC 
(and also for ECM=7 TeV as a reference). In our case, the cross-section for the light 
(heavy) Higgs boson will be scaled by cos 2 a (sin 2 a) compared to the total Higgs pro- 
duction cross-section which is shown in Fig. |3|. Fig. || shows the variation of the decay 
branching fraction of the right-handed neutrino, ^f, with its mass for the fixed values of 
Higgs masses (m^ = 125 GeV, m# = 300 GeV), cos a = 0.8 and A couplings chosen for 
the case of BP4, which is defined later. We can read from the plot that the for low mass of 
the right-handed neutrino, it mostly decays to gauge bosons but for mq, > 250 the decay 
branching to hv dominates (always > 50%). 

Let us now select some relevant points for the phenomenological studies. For this 
purpose, we choose sin 9 sa mp/Mpi = 0.05 close to the current bound ( jO| ) for each 
generation. With this choice, we further consider two different options for the light Higgs 
boson. 

1. Very low mass light Higgs (h): A very light (non-standard) Higgs boson mass (below 
114 GeV) can still satisfy the corresponding LEP bound as discussed in the previous 
section. Keeping the the light Higgs mass fixed at 50 GeV we now choose two sets 
of points: 

(m H ,mq,,c a ) = (125 GeV, 100 GeV, 0.1/0.25). 

Thus, for these two points, it is the heavier Higgs which stays in the discovery region 
of the LHC. 

2. Light Higgs (h) in the LHC discovery region: Here we take mh = 125 GeV and choose 
the following combinations: 

(m H ,m*,c a ) = (200/300 GeV, 100 GeV, 0.8). 

Note that the LHC bound (Fig. |l|) sets a limit of s 2 < 0.33(0.36) for m H = 200(300) 
GeV. So, we made a generous choice of c a = 0.8 corresponding to s 2 = 0.36. 

Given the Higgs masses and (v, v') = (246, 12000) GeV, the four-point Higgs couplings 
for each benchmark point in Table |] can be found as follows: 

1. BP1: Ai = 0.128, A 2 = 9.14 x 10" 6 , A 3 = 4.42 x 10~ 4 , 
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Benchmark 
Points 


mfj 
(GeV) 


rriH 
(GeV) 


(GeV) 


cos a 


BP1 


50 


125 


100 


0.1 


BP2 


50 


125 


100 


0.25 


BP3 


125 


200 


100 


0.8 


BP4 


125 


300 


100 


0.8 



Table 1: Benchmark points for common value of sm8 = 0.05. 



2. BP2: Ai = 0.122, A 2 = 1.15 x 10~ 5 , A 3 = 1.08 x 10~ 3 , 

3. BP3: Ai = 0.202, A 2 = 1.08 x 10~ 4 , A 3 = 3.96 x 10~ 3 , 

4. BP4: Ai = 0.350, A 2 = 2.20 x 10" 4 , A 3 =1.21 x 10~ 2 . 

Table | summarizes the benchmark points for our analysis. From Table |], we can see that 



BP1 and BP2 correspond to low mass of the light Higgs that is not excluded by LEP [ 13 [ 
and the heavier Higgs mass is in the discovery region of the LHC. BP3 and BP4 correspond 
to relatively heavier light Higgs = 125 GeV and heavy Higgses are 200 and 300 GeV, 
respectively. 

In Table g[ the decay branching fractions of the heavy Higgs boson H are shown for 
the benchmark points. For BP1 and BP2, the heavy (125 GeV) Higgs decays to is not 
open because of the unavailable phase space, but the decays to the light Higgs boson (h) 
pair is open. In case of BP3 where mn = 200 GeV, the H decays to gauge bosons via one 
off-shell gauge boson as before, but the decay modes to and hh are closed. In case of 
BP4, the heavy Higgs boson can decay to \t"I', as well as to hh. Table |3| shows the decay 
branching fractions for the light Higgs boson h for all the benchmark points. While the h 
decay to is closed in case of BP1 and BP2, it is open for BP3 and BP4 with a branching 
fraction ~ 7%. Table || gives the decay branching fraction of for the four benchmark 
points. From the Table ||, it is clear that the right-handed neutrino mostly decays to gauge 
bosons. 

In Table ^, we present the cross-sections of the Higgs boson production for the bench- 
mark points with the center of mass energy of 14 and 8 TeV at the LHC, where we used 
CTEQ5L as PDF and v^asa scale. 
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Decay Modes 


Branching Fraction 




BP1 
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U.UUZ 




n nf; 

U.UD 


U.Uo 


n nni o 

U.UU1Z 


y.i x iu 


bb 


n K7 
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U.Z 1 


n nn/i 

U.UU4 
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U.UD 


n m 
U.Uo 


4.U X 1U 


9 i v i n — 4 
Z.l X iu 


hh 


0.035 


0.55 




0.63 


WW/WW* 


0.19 


0.09 


0.629 


0.20 


zz/zz* 


0.017 


0.009 


0.364 


0.16 


gg 


0.046 


0.021 


0.001 


0.002 


77 


0.002 


0.001 


1.5 x 10" 4 


7.6 x 10" 5 



Table 2: Decay branching fraction of H for four benchmark points. 



Let us now we discuss the various hnal sates for further phenomenological studies of 
the benchmark points. We focus on the contribution of the right-handed neutrino from the 
decays of heavy and light Higgses depending on the benchmark points. The final states 
are then determined from the decays of the right-handed neutrino. 

4.1 BP1 & BP2 

These two benchmark points correspond to the case of nih = 50 GeV, forbidding the decay 
h — > v^. The main decay modes of the lighter Higgs boson are h —> bb and rr leading 
to the final states: 2-6-jets + no $t and 2-r-jets + no As the case with the usual 
Higgs search, it is very difficult to isolate the former signal from the QCD backgrounds. 
On the other hand, the r-jet analysis could be interesting. As the lighter Higgs is very 
light in our case, it leads to very soft jets and thus the r-jet tagging efficiency goes down. 
However, r-jet tagging through one prong decay of r can have better handle over the QCD 
background. Thus, a possibility of 2-r-jets + no ]/t final state can still be important with 
a large accumulated luminosity. 

The non-standard Higgs signature of our interest comes from the heavy Higgs which is 
SM-like with rriH = 125 GeV. As can be seen from Table ^ with mq, = 100 GeV, the decay 
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Decay Modes 


Branching Fraction 




BP1 


BP2 


BP3 


BP4 








U.U f O 


U.U 1 O 


bb 


U.o 1 


U.o 1 DO 


U.Oo 


U.Oo 


TT 


0.076 


0.076 


0.06 


0.06 


WW 






0.19 


0.19 


zz 






0.017 


0.017 


gg 


0.011 


0.011 


0.047 


0.047 


77 


0.003 


0.003 


0.002 


0.002 



Table 3: Decay branching fraction of h for four benchmark points. 



Decay Modes 


Branching Fraction 




BP1 


BP2 


BP3 


BP4 


hu 


0.068 


0.167 






Zu e 


0.116 


0.103 


0.124 


0.124 


W + e 


0.817 


0.730 


0.876 


0.876 



Table 4: Decay branching fraction of for four benchmark points. 

H — > vty has a branching fraction ~ 6% for BP1 and 3% for BP2. Then, the produced 
right-handed neutrino (^f) will decay to gauge bosons as we can read from Table |j. This 
will lead to the prominent signal with respect to the backgrounds (see Fig. |5|): 

21 + 0-jet+ fa. (4.1) 

Hadronically quiet criteria along with two leptons kills the dangerous QCD and other 
hadronic backgrounds. The major SM hadronically quiet backgrounds come from the 
leptonic decays of WW, WZ and ZZ. Although there is a huge single vector boson 
production leading to two leptons; pp — > "y*/Z — > £ + £~, this does not carry any missing 
energy and can be efficiently reducible, specially with the Z mass window cut as will be 
discussed later. 
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ECM=14 TeV 


ECM=8 TeV 




0* (fb) 


a H (fb) 


<? h (fb) 


a H (fb) 


BP1 


872.33 


15869.3 


370.18 


5601.02 


BP2 


5452.06 


15027.75 


2313.63 


5304.0 


BP3 


10258.94 


2515.64 


3620.86 


782.89 


BP4 


10258.94 


1386.12 


3620.86 


376.92 



Table 5: gg —> h/H production cross-section for four benchmark points for ECM=14 and 8 TeV, 
respectively, using CTEQ5L as PDF and \f§ as scale. 




Figure 5: Decay topology of the Higgs(es) leading to the opposite sign lepton pair in the final 
state. 

Let us remark that, in the case of BP2 with cos a = 0.25, the heavy Higgs (H) decay 
branching fraction to light Higgs (h) pair is about 55%. This will generate a 4b final state 
with no pr- Plotting the invariant mass distribution of 6-jet pairs, one could find a h peak 
although the signals would be challenged by the QCD background. 

4.2 BP3 & BP4 

For both BP3 and BP4, the light Higgs boson {h) can decay to vty with branching fraction 
7.3%. This leads to a hadronically quiet di-leptonic final state (|4.1| ), which was the case 
with the heavy Higgs for BP1 and BP2. The hadronically quiet 11 scenario thus a generic 
signature for all the benchmark points. 

On the other hand, the heavy Higgs {H) decay branching fraction to v^> drops to 
< 1%. In case of BP4, the heavy Higgs decay to opens up but with very small 
branching fraction ~ 0.2%. Even with such small branching fractions, the and Vfy 
modes can still be probable due to the large production cross-section as in Table ||. BP4 
allows the heavier Higgs decay mode H — > hh with 63% of decay branching fraction, which 
is closed for BP3. The subsequent light Higgs decay to will open up an additional final 
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state: 

U + O-jet + fa. (4.2) 

The At final state in case of benchmark points 3 and 4 also comes from the heavier Higgs 
decay to gauge boson pair. This implies that for heavier Higgs scenarios, either through 
gauge boson pairs or through the right-handed neutrino pair, 2£/4£ final states without 
any hadronic activity are generic ones. In the following section, we will focus on these 
hadronically quiet two lepton final states for a collider simulation as the four lepton states 
turn out to be less significant. 



5. Results 

In this section, we go through a PYTHIA level analysis before presenting the final state 



results. PYTHIA (version 6.4.22) [ ]15p has been used for the purpose of event genera- 
tion. The model is implemented in Calchep |]l6| and the corresponding mass spectrum and 
decay branching fractions are generated. These are then fed to PYTHIA by using the SLHA 
interface [flTfl . Subsequent decays of the produced particles, hadronization and the collider 
analysis were performed using PYTHIA. We used CTEQ5L parton distribution function (PDF) 
|l8tl for the analysis. The renormalization/factorization scale Q was chosen to be the par- 
ton level center of mass energy, y/~§. We also kept ISR, FSR and multiple interaction on 
for the analysis. We have used PYCELL, the toy calorimeter simulation provided in PYTHIA, 
with the following criteria: 

I. The calorimeter coverage is \rj\ < 4.5 and the segmentation is given by Arj x A<j> = 

0. 09 x 0.09 which resembles a generic LHC detector. 

II. AR = ^{Ar/) 2 + (A(f>) 2 = 0.5 has been used in cone algorithm for jet finding. 
HI 4%n = 20 GeV. 

IV. No jet matches with a hard lepton in the event. 

In addition, the following set of standard kinematic cuts were incorporated throughout: 

1. p e T >5 GeV and \q\i < 2.5, 

2. \rj\j < 2.5, AR £j > 0.4, AR U > 0.2, 

where ARgj and ARn measure the lepton-jet and lepton-lepton isolation, respectively. 
Events with isolated leptons, having px > 5 GeV, are taken for the final state analysis. 

Fig. |6] (left) plot shows the j/t distribution for the model for BP4 and the SM back- 
ground which is coming from the gauge boson pairs (WW, ZZ, ZW). In both cases, the 
origin of $/t is neutrinos. In the case of background, it comes from the decay of the gauge 
bosons and thus it peaks around 45 GeV. On the other hand, in the case of signal, there are 
two different sources of neutrinos. The ones from the gauge bosons will behave similarly 
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Figure 6: $t distribution (left) from heavy (H) and light (h) Higgses for BP4 and gauge boson 
pair background. The lepton px distribution (right) of signalxl20 for BP4 heavier Higgs (H), 
signal x 10 for the lighter Higgs (h) and background at 10 fb _1 integrated luminosity. Both of the 
plots are for ECM=8 TeV. 

but the ones from the Higgs decay (H — > ^Sf v) can have more pr, depending on the mass 
difference between the Higgs and the right-handed neutrino. This results in a tail at higher 
j/t ( see Fig. H (left)) and some enhancement in the j/t distribution for the heavier Higgs 
signal around 100 GeV for the case of benchmark point 4. In comparison, for the lighter 
Higgs (h), the tail dies below 100 GeV as expected. A requirement of minimum jfr could 
be a good handle to kill the leptons coming from single Z, as charged lepton pair and 
neutrino pair come from Z decay in mutually exclusive manner. Higher pr cut > 100 GeV 
could be a good handle to search for the existing heavy Higgs (H), however, we use soft 
pr cut while looking for 125 GeV Higgs as discussed in the next section. 

Fig. U (right) plot shows the lepton pt distribution coming from the signal (H and h) 
in the case of BP4 and from the gauge boson pair background, respectively. In the case 
of signal, the leptons that we will be looking for, have different origin. The one coming 
from the right-handed neutrino can be very soft due to the small mass gap between the 
right-handed neutrino and the W boson. The other one coming from W decay would be 
of SM like. We can see from the figure that the lepton from the signal, specially for the 
case of the heavier Higgs (H), is as hard as 125-150 GeV. This is because for BP4 the 
heavier Higgs mass is 300 GeV. The decay of such a particle to the light Higgs (h) pair 
or gauge boson pairs or though tyty or will share the momentum. Now the boosted 
decay products will transfer their momentum to their daughter leptons, which carry the 
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~~ ! ' ' ' BP4h b : 

■ VV I 

1 'VV ^ 

E 

1 2 3 4 5 6 7 
Number of leptons e-|x 

Figure 7: The left figure shows the lepton multiplicity distribution of signal for BP4 heavier Higgs 
(H), for the lighter Higgs (h) and background at 10 fb _1 integrated luminosity with ECM=8 TeV 
and the right figure shows n e — (n^ — n e ) from the signal and di-boson background denoted by 
VV ('VV'). 

momentum depending on the decay channels. On the other hand, the lepton coming from 
the light Higgs (h) decay could be of very small momentum due to the small mass difference 
between the Higgs and the right-handed neutrino, which decays further to a lepton and 
W. Thus, the lepton pr cuts, upper for the soft leptons and lower for hard leptons, will be 
crucial for signal event selection and reduction of the backgrounds. 

Fig. |7| (left) describes the lepton multiplicity distribution for the signal for BP4 com- 
ing from the heavier Higgs and the lighter Higgs and also from the gauge boson pair 
backgrounds, respectively. The lepton multiplicity distribution is independent of the rel- 
ative charge, i.e., same sign or opposite sign. We can see from the Fig. ^ (left) that the 
gauge boson pair backgrounds also have 21 and 4/ final states from their decay to the 
charged leptons. Fig. |7| (right) describes the electron and muon number difference. The 
signal events contain more electrons for the signal from the decay of the right-handed neu- 
trino as it is assumed to couple only to one flavour, i.e, the electron. On the other hand, 
for the background we expect the number to be very similar for e and /x case (that is, 
\n e — n^\ <C n e + n^). From Fig. ^ (left and right), we can see that the background is 
reduced a lot, in particular, for the case of n e — n^. 

The main Standard Model backgrounds come from the di-boson production as dis- 
cussed earlier. Apart from these, the Drell-Yan, i.e., pp — > Z/'y* — > £ + l~ could have 
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SI 
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h 
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0.09 
352.02 


0.23 
154.34 


269.03 
43.75 


269.03 
4.60 


3558.47 


576.89 


209.29 


Significance 




5.1 
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(20 


h 
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42.88 
820.44 


267.97 
324.60 


600.66 
117.48 


600.66 
14.34 


5774.59 


3388.27 


546.85 


Significance 




8.4 


5.8 


7.0 


6.0 









Table 6: Number of events for 2£ final states for the benchmark points and the SM backgrounds 
at an integrated luminosity of 10 fb _1 with ECM = 8 and 14 TeV. 



been a major background due to large cross-section but we find that this background is 
reducible one as it gives dideptons without missing energy. We will also have a window 
cut \Mit — Mz\ > 5 GeV to further kill the dominant Z induced background. The ti 
background which is sub-dominant but still comparable to the signal event has also been 
implemented. 

In the following two subsections, we present two sets of results. First, we discuss 
the hadronically quiet di-lepton scenario for inclusive flavour. Then, we address the case 
of lepton flavour violation, that is, y e 3> y^ ]T (or y^ S> y e ,r) for which the right-handed 
neutrino will decay only to eW and v e Z (or to [iW and v^Z). 

5.1 Inclusive charged lepton signature 

In this section, we discuss the flavour independent results which include both electron 
and muon, though we have lepton flavour violating coupling as the right-handed neutrino 
couples mostly to electrons, i.e., y e 3> y^r- This study is also suitable for the case of 
Ve = y^ 

For BP1 and BP2, h — > decay is closed and h dominantly decays to bb or tt (see 
Table |3|). The only contribution to 21 comes from the semi-leptonic decays of b and r. 
Thus, we do not have large contribution from the light Higgs in these cases. BP2 has 
more events than BP1 in the final state due to the larger cross-section compared to BP1. 
On the other hand, the heavier Higgs (H) can decay to with small branching fraction 
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(< 10%) as described in Table Thus, the contributions are large enough due to large 
production cross-sections for the heavier Higgs in the cases of BP1 and BP2 as can be seen 
from Table ||. However, in the final event counting we expect BP1 should have twice the 
event of BP2, as the Bi(H — > ^u) ~ 6% in the case of BP1 and which is 3% in BP2 as 
cane be read from Table 

BP3 and BP4 are exactly same in terms of the lighter Higgs (h) having = 125 GeV 
and c a = 0.8. This leads to the same decay branching to i/^f, 7.3%, which results in large 
number of events for the final state. Also the gauge boson pairs contribute to the final 
states through off-shell decays of the lighter Higgs (h). In the case of the heavier Higgs (H), 
BP3 has the cross-section almost twice as large as BP4 (see Table |||), and the branching 
fraction to gauge boson pair is around ~ 99% for BP3 and ~ 36% for BP4. Unlike BP3, 
BP4 allows the decay H — > hh with the branching fraction of 63%, which gives additional 
contribution to the di-lepton final state. 

As we have discussed earlier, the lepton coming from the right-handed neutrino will 
be softer compared with that coming from (see Fig. |6|). Thus, we demand the softer 
lepton with pt < 30 GeV and the harder lepton with pj- > 20 GeV. This would help to 
reduce the SM di-lepton backgrounds, which generically come from gauge-boson decays. 
To kill the Z boson background , we reject the events with lepton invariant mass (M^) 
around Z mass window, i.e., — Mz\ < 5 GeV (defined as nz = 0). On top of that, we 
also demand p^ + p^p < 100 GeV and pr > 30 GeV. When we look for the signal coming 
from 125 GeV Higgs, we demand the sum of lepton px and $t should be less than 125 
GeV. This cut will exclusively select events coming from 125 GeV Higgs. When we look 
for the events coming directly from heavier Higgs (H: mn = 200,300 GeV), the cut can 
be modified accordingly. In summary, the selection cut used for the inclusive 21 analysis is 
as follows: 

SI : n e = 2, n jets = 0, n z = 0, r/ T > 30GeV, P^ < lOOGeV, M eff < 125GeV (5.1) 

where M e g- = YKPt~^~ Pt) an d fj> = YIPt- 

In Table ^ we present all the event numbers from the signal (S) and background (B), 
and the significance defined by S/y/ S + B for an integrated luminosity of 10 fb _1 at the 8 
and 14 TeV LHC. We can see that the behaviour of BP1 and BP2 are complementary to 
that of BP3 and BP4 as the main contribution to the signal events come either from the 
heavy Higgs (H) for BP1 and BP2 or from the light Higgs (/i) for BP3 and BP4, and the 
sum of two contributions are more or less the same except for BP2. The signal significance 
at the 8 TeV LHC reaches around 5a except for BP2. In the case of the 14 TeV LHC, all 
the benchmark points reach high significance at an integrated luminosity of 10 fb , which 
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0.04 
122.38 


0.0 
55.96 


106.82 
10.44 


106.82 
1.58 


878.24 


263.3 


59.03 


Significance 




3.4 


1.6 


3.2 


3.0 
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BP4 


VV 




tt 


S2 

(2e) 


h 
H 


11.00 
311.83 


68.70 
125.48 


242.62 
26.67 


242.62 
3.81 


1480.36 


1602.50 


123.78 


Significance 




5.4 


3.3 


4.6 


4.2 









Table 7: Number of events for 2e final states for the benchmark points and the SM backgrounds 
at an integrated luminosity of 10 ftr 1 with ECM=8 TeV and 14 TeV. 



implies that the parameter space of the inverse seesaw model can be readily probed at the 
LHC. 

5.2 Lepton Flavour Violating signature 

As mentioned earlier, much more significant search can be made if the right-handed neu- 
trino has hierarchical Yukawa couplings and thus allows lepton flavour violation in its decay. 
Let us first consider the case that the right-handed neutrino couples only to electron. A 
comparative study of electron and muon lepton flavour will give a vital clue about the 
structure of the model as it can reduce significantly the SM background. 

For this study, we select each flavour state keeping all the other cuts the same as in the 
previous section (^^), and present results of our analysis for each signal channel labeled as 

S2:2e; S3 : 2/i ; S4 : le + lfi ; and S5:2e-2/^. (5.2) 

As expected, S5, corresponding to S2-S3, will give a novel signature to probe the inverse 
seesaw mechanism. Instead of S5, one could make an equivalent study with other lepton 
flavour violating final states such as 2 • S2 — S4 or 2 • S3 — S4. 

Table |7| (||) shows the number of electron (muon) events for the signal for all the 
benchmark points and the backgrounds at an integrated luminosity of 10 fb" 1 with center 
of mass energy of 8 and 14 TeV at the LHC. Unlike the lepton-universal gauge boson 
decay, the right-handed neutrino will produce only electron in the IW mode. For BP1 and 
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894.30 


290.48 
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Significance 




1.4 


0.65 


1.1 
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Significance 




2.0 


2.0 


1.5 


1.0 









Table 8: Number of events for 2/j, final states for the benchmark points and the SM backgrounds 
at an integrated luminosity of 10 fhr 1 with ECM=8 and 14 TeV. 



BP2, h — > decay is not open and thus we do not have the extra contribution to the 
electron events in the case of the light Higgs boson (h). This is clear from Table [7| and 
Table |8| where the electron and muon event numbers are similar as they come from the 
semi-leptonic decays of the decay products of the light Higgs (h). On the other hand, the 
heavy Higgs (H) with run = 125 GeV, can decay to vty and thus produces more electrons. 
The BP1 contribution is more than BP2 due to the larger branching fraction to Vfy and 
the cross-section of the former. 

Similar is the case with the lighter Higgs (h) contribution for BP3 and BP4 as estab- 
lished in Table |?] and |8| The situation for BP3 is again a little different, as the heavier 
Higgs decay H — )• hh is not allowed and H — > z/\T/ branching fraction is very small ~ 0.1% 
(see Table |2|). This ceases the extra electron contribution over muon for BP3. For BP4, 
heavier Higgs mass is 300 GeV and thus H —> hh is allowed with 63% branching fraction. 
The lighter Higgs thus produced will decay to the right-handed neutrino as can be read 
from Table ||. Similar to the previous 11 case, the complementary behaviour between the 
heavy and the light Higgs remains, i.e., for BP1 and BP2, it is the heavier Higgs which gives 
the lepton flavour violating signal generating extra electrons in the final state, whereas for 
BP3 and BP4, it is mainly the lighter Higgs. From Table |?] and ||, we can see that the 
significance for S3 is much smaller compared to S4 as expected. 

For completeness, we also present the event numbers for the final state with le + 1/i, 
defined as S4. Table || shows the event numbers for S4. The Drell-Yan background is 
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Table 9: Number of events for le+l/i final states for the benchmark points and the SM backgrounds 
at an integrated luminosity of 10 fb" 1 with ECM=8 and 14 TeV. 
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Table 10: Number of events for 2e — 2/i final states for the benchmark points and the SM back- 
grounds at an integrated luminosity of 10 fb -1 with ECM=8 and 14 TeV. 



reduced a lot due to the demand of different flavour leptons in the final state. The signal 
significance in this case is comparable to S2 (2e) and greater than S3 (2fi). 

S2 (2e) from Table [7| and S3 (2[i) from Table [8| lead us to check the difference in the 
event numbers for the 2e and 1\i final states. As the SM backgrounds come from the 
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Figure 8: The opposite sign lepton invariant mass distribution with events with the same cut as 
Sll for BP1 heavier Higgs (if) at 10 ftr 1 integrated luminosity with ECM=8 TeV. 



flavour-blind decays of the gauge bosons, we expect that 2e — 2/i event number will kill 
the background substantially giving only the extra electrons coming from the right-handed 
neutrino decays. This is the artifact of the lepton flavour violating right-handed neutrino 
coupling: y e 3> y^, T - Table [H] shows the difference in number of 2e and 2[i events (S5). 
We can again see that 125 GeV Higgs, which is the heavy Higgs (H) for BP1 and BP2 
and the light Higgs (h) for BP3 and BP4, decays into more electrons than muons. The 
negative sign shows the events with more muons than electrons. From Table 10, we can see 
that 2e — 2/i signal has a high significance except for BP2 for which the signal events are 
overshadowed by background for both 8TeV and 14TeV LHC at an integrated luminosity 
of 10 fb _1 . 

Based on the above results, we can get the results for the opposite case where the 
right-handed neutrino couples only to the muon flavour: y^ 3> y e ^ T . In this case, the signal 
significance drops due to the muon excess in the background: 7 a (Wa) for all points except 
BP2 with the 8 (14) TeV LHC which is still better than the inclusive flavour searches. 



5.3 Mass measurement 

Let us make a comment on the prospect of measuring the right-handed neutrino mass 
through the di-leptonic edge. Since we are looking for the decay chain of the right-handed 
neutrino; Vl/ — > IW — > 11' v' , the invariant mass of the final leptons gives rise to the famous 
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di-lepton edge [19]: 



mTr = m n l-^Jl-^. (5.3) 

As m„ rs 0, the measurement of the maximum di-lepton invariant mass will tell us about 
the value of my. 

Fig. H shows the distribution of the di-lepton invariant mass for opposite sign di-leptons 
from the signals only under the same cuts as used for SI. Clearly, the edge can be seen at 60 
GeV. Now reconstructing W through its hadron decay mode would be crucial in detecting 
the decay topology. For the signal, the decay topology via W boson indicates the mass of 
the right-handed neutrino at 100 GeV from the Eq. (|5.3|) . Thus, the prescription described 
above not only gives a variable to get signal events over backgrounds but also results in 
a possibility to measure the mass of the right-handed neutrino. Analyzing the signal and 
background for events at 10 fb _1 of integrated luminosity, however, we find that it is very 
difficult to recognize the di-leptonic edge. We have to go to much higher luminosity to get 
a clear edge from the total distribution. 



6. Conclusion 

The inverse seesaw model introduces a tiny B — L breaking Major ana mass of the right- 
handed neutrino which explains the smallness of the neutrino mass. This allows rather large 
neutrino Yukawa couplings through which the right-handed neutrinos can be produced at 
the LHC. We point out that the di-lepton final state with missing energy could be a smoking 
gun signal probing the Higgs and the right-handed neutrino of the inverse seesaw model. 
Furthermore, a novel signature of lepton flavor violation in these final states, such as the 
difference in the ee and /i/x event numbers, may occur due to flavor-dependent neutrino 
Yukawa couplings. Taking the neutrino Yukawa coupling y e ^ = 0.029, which is close to 
the upper limit put by the lepton universality, and the right-handed neutrino mass of 
100 GeV, we studied the LHC prospects to look for the inverse seesaw model in the four 
benchmark points. Performing a PYTHIA level simulation, it is found that the 5a signal 
significance can be achieved with the integrated luminosity of 10 — 20 fb _1 for the inclusive 
lepton (flavour-blind) signature, and ^ 2 fb _1 for the flavour violating signature at the 8 
TeV LHC. We also pointed out that the observation of the di-lepton edge could shed light 
on the right-handed neutrino mass measurement. But it turns out to be hard to see the 
di-leptonic edge over the background with the nominal luminosity of the LHC14. 
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A. Vertices 



We write down all the relevant vertices below: 



n ~ Qu(d), ~ Quid), ■ ^= Pr 



(A.l) 



H - q u ( d ) t ~ Quid), ■ ™ (d> ^2 P R> ( A - 2 ) 

h-k-i,: V -^Pn, (A.3) 

H-U-U: yj ^Pn> (A.4) 

h-^-n : ^osacos^sinasin^ ^ 

P - % - n : ^ SinaC ° S ^ COSaSin V L , (A.6) 

- T yw. cos a sin 6* — y s . sin a cos „ , , _. 

: ^ Pl, (A.7) 

~ T y u . sin a sin 6* + y s . cos a cos „ . 

: ^ ^ P L , (A.8) 

v2 

WH--^-^ : ^^P L , (A.9) 

: ^^k p L , (A.10) 
V2 

*-*-».: f^ft, (A.11) 
2 cos 0vf 

Z : f^P L , (A.12) 
2 cos 6>vk 

„ T / „ t <7 cos sin „ . . . 

Z-Pi-^ =^-^-1/0 : ^ P ^ ( A - 13 ) 

2 cos 
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h 4 : 6[\ 1 c 4 a + X 2 s 4 a + X 3 c 2 a s 2 a ], (A.14) 

H 4 : 6[Ai4 + A 2 4 + A 3 4*a], ( A -!5) 

h 2 H 2 : 6(Ai + M)c 2 a s 2 a + A 3 (4 + <£ " 4 44), (A.16) 

/i 3 tf : (6A1 - 3As)4aa + (-6A 2 + 3A 3 )Ca4. (A. 17) 

/1F 3 : (-6A2 + 3A 3 )c 3 s Q + (6A1 - 3A 3 )c a c 3 , (A.18) 

/i 3 : 6Ait;c 3 — 6A2f's 3 + 3A 3 (t;s Q , — f'c a )s a c Q , (A. 19) 

F 3 : 6Aiws 3 +6A2WC 3 +3A 3 (uc Q + v's a )s Q c Q , (A. 20) 
hH 2 : 6Xivc a s 2 t -6\2v's a c 2 l + \3[2(-vs a + v'c a )s a c a -v's^ + vc^], (A. 21) 
/i 2 F : Q\ivc 2 a s a + §\ 2 v's 2 a c a + \z[-2(v's a + vc a )s a c a + v'cl + vsll (A.22) 

WW/} : 2m 2 v c a /v, (A.23) 

WWH : 2m^s Q /?;, (A.24) 

WWhh : 2m^c^/u 2 , (A.25) 

WWhH : 2m^s Q c Q /t; 2 , (A.26) 

WWHH : 2m 2 v s 2 a /v 2 , (A.27) 

ZZ/i : 2m 2 z c a /v, (A. 28) 

: 2m|s Q /t>, (A. 29) 

ZZhh : 2m 2 z c 2 a /v 2 , (A.30) 

ZZhH : 2m 2 z s a c a /v 2 , (A. 31) 

: 2m 2 z s 2 a /v 2 , (A.32) 

Z'Z'fc : -2|y*_J 2 m|,^, (A.33) 

Z'Z'tf : 2|Y£_ L | 2 m|,^, (A.34) 

Z'Z'/i/i : 2|y£_ L | 2 m|,^§, (A.35) 

Z'Z'hH : -2\Y^ L \ 2 m 2 z ,s a ^, (A.36) 

Z'Z'HH : 2\Y*_ L \ 2 m 2 z ,^, (A.37) 
where Y£_ L = —1/2 in our case. 

References 



[1] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716 (2012) 30 [arXiv:1207.7235 
[hep-ex]]. 



- 24 - 



[2] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716 (2012) 1 [arXiv:1207.7214 [hep-ex]]. 

[3] R. N. Mohapatra, Phys. Rev. Lett. 56, 561 (1986). 

[4] R. N. Mohapatra and J. W. F. Valle, Phys. Rev. D 34, 1642 (1986). 

[5] S. Khalil, J. Phys. G 35, 055001 (2008) [arXiv:hep-ph/0611205]. 

[6] W. Abdallah, A. Awad, S. Khalil and H. Okada, Eur. Phys. J. C , 72:2108 (2012) 
[arXiv:1105.1047 [hep-ph]]. 

[7] P. S. B. Dev, R. Franceschini and R. N. Mohapatra, arXiv:1207.2756 [hep-ph]. 

[8] C. G. Cely, A. Ibarra, E. Molinaro and S. T. Petcov, arXiv:1208.3654 [hep-ph]. 

[9] CMS Collaboration, Phys. Lett. B 710 (2012) 91 [arXiv:1202.1489 [hep-ex]]; ATLAS 
Collaboration, arXiv:1206.0756, ATLAS-CONF-2012-098. 

[10] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[11] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 

[12] F. del Aguila, J. dc Bias and M. Perez- Victoria, Phys. Rev. D 78, 013010 (2008) 
[arXiv:0803.4008 [hep-ph]]. 

[13] R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH Collaboration 
and and], Phys. Lett. B 565, 61 (2003) [arXiv:hep-ex/0306033]. 

[14] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710 (2012) 26 [arXiv:1202.1488 
[hep-ex]]. 

[15] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 026 (2006). 

[16] A. Pukhov, hcp-ph/0412191. 

[17] P. Z. Skands et al., JHEP 0407, 036 (2004). 

[18] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000); J. Pumplin et al, 
JHEP 0207, 012 (2002). 

[19] C. G. Lester and D. J. Summers, Phys. Lett. B 463 (1999) 99 [hep-ph/9906349]. 



- 25 - 



